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A ROBUST MULTIGRID METHOD FOR THE TIME-DEPENDENT 

STOKES PROBLEM* 

STEFAN TAKACSt 

Abstract. We propose a coupled multigrid method for generalized Stokes flow problems. Such 
problems occur as subproblems in implicit time-stepping approaches for time-dependent Stokes prob¬ 
lems. The discretized Stokes system is a large-scale linear system whose condition number depends 
on the grid size of the spatial discretization and of the length of the time step. Recently, for this 
problem a coupled multigrid method has been proposed, where in each smoothing step a Poisson 
problem has to be solved (approximately) for the pressure field. In the present paper, we propose 
a coupled multigrid method where the solution of such sub-problems is not needed. We prove that 
the proposed method shows robust convergence behavior in the grid size of the spatial discretization 
and of the length of the time-step. 
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1. Introduction. We consider the following model problem (generalized Stokes 
flow problem). Let C be a bounded polygonal domain and assume / S 
and g G L^(fi) to be given. Find a velocity field u and a pressure distribution p such 
that 


—Au + /3u + Vp = / in n, 
V ■ u = g in ft, 
rt = 0 on dfi 


( 1 . 1 ) 


is satisfied. /I > 0 is assumed to be a given parameter. The problem (1.11 appears as 
an auxiliary problem for implicit time-stepping approaches to solve an incompressible, 
time-dependent Stokes flow problem. In this case, the parameter (3 is proportional to 
the inverse of the length of the time-step, scaled by a viscosity parameter. 

To obtain existence and uniqueness of the solution, we further require f^pdx = 
J^gdx = 0. Note that the analysis presented in this paper is (due to the use of 
regularity results) restricted to convex domains. However, the numerical method 
presented in the paper can be applied also to non-convex domains. 

The discretization of the problem leads to an indefinite linear system with saddle- 
point structure. The main goal of this work is to construct and to analyze numerical 
methods that produce an approximate solution to the problem, where the computa¬ 
tional complexity can be bounded by the number of unknowns times a constant which 
is independent of the grid level (of the spatial discretization) and the choice of /3, in 
particular for large values of (3 (which correspond to small time steps). 

For the solution of such a saddle-point problem, there are several possibilities. 
In 0 [TTJ [T31HH [TSJ [5^ , various kinds of preconditioners have been proposed for this 
problem which can be combined with a Krylov-subspace method as outer iteration 
scheme to yield an iterative solver for the problem. 

An alternative is to apply a multigrid algorithm directly to the coupled system. 
On the one hand such methods are typically quite fast, on the other hand such an 
approach does not require an outer iteration scheme. For /3 = 0, problem (1.11 is 
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the standard (stationary) Stokes problem. For this case several multigrid solvers are 
available, see, e.g., nzuniiiaiMi and the papers cited in Ellis]. The construction 
of a multigrid method for /? > 0, particularly if the method is desired to show robust 
convergence behavior in /3, is more involved, see m for an overview and numerical 
results. Recently, an important step forward has been archived by Olshanskii, who 
has proposed such a robust multigrid method, see |15| . In the named paper however, 
a Poisson problem for the pressure has to be solved approximately (for example by 
applying one V-cycle, cf. Section 8 in m) for each step of the smoothing iteration. 
(We will comment on this in Remarks |4.1| and 4.41. 

The goal of the present paper is to drop this requirement and, therefore, reduce 
the computational costs. We propose a multigrid solver where the smoother is a simple 
linear iteration scheme and we prove that the proposed method is robust in the grid 
size of the spatial discretization and in the choice of /3. We will present a convergence 
proof for our multigrid method based on the classical splitting of the analysis into 
smoothing property and approximation property, see E|- 

The results of this paper form a basis of the convergence analysis of an all-at-once 
multigrid method for the Stokes optimal control problem, cf. |19| . 

This paper is organized as follows. In Section]^ we will introduce the variational 
formulation and discuss its discretization. The multigrid framework will be presented 
in Section]^ In Sectionj^we will discuss the choice of the smoother. The proof of the 
approximation property will be given in Section]^ Numerical results which illustrate 
the convergence result will be presented in Section]^ In Section]^ we will close with 
conclusions. 

2. Variational formulation and discretization. Here and in what follows, 
L^(n) and denote the standard Lebesgue and Sobolev spaces with associated 

standard norms || • ||l 2 (q) and || • H^qn), respectively. The space To(^) the space 
of all functions in Lp‘{Vl) with mean value 0. The space Hq{^) is the space of all 
functions in H^{^) that vanish on the boundary. Both spaces are equipped with 
standard norms, i.e., || • ||i2(f2) := || • 11^2(0) and || • |lffi(n) := II • llggo)- 

Using these spaces, we can set up the variational formulation of (1.11, which reads 
as follows. Find u G U := [FIo(f2)]^ and p G P := To(^) such that 


(Vu, Vu)i2(n) + /?(«, u)l2(q) + (p, V • M)L2(a) = (/, m)l2(q) 

(V • u,p)l2{q) = {g,p)mQ) 

holds for all u G U and p G P. Certainly, the variational problem can be rewritten as 
one variational equation as follows. Find x G X such that 


B{x, x) = P{x) for all x G X, 


( 2 . 1 ) 


where X := U x P, x := {u,p), x := {u,p) and 

B{{u,p), iu,p)) := (Vu, Vu)i 2 (Q) +/3(u,{i)i2(Q) + (p, V • u)L^(n) + (V • u,p)i 2 (Q), 
P{u,p) := if, u)L^(n) + (p,p)L2(n). 


We are interested in finding an approximative solution for equation (2.1). The 


convergence analysis follows standard approaches, i.e., we show that the problem in 
question is well posed in some norm || • ||x (which is for the Poisson problem the H^- 
norm and for the standard Stokes problem the pair of for the velocity and for 
the pressure). In our case this norm will be parameter-dependent. In a second step. 
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we will introduce a properly scaled L^-like norm ||| • |||fe, which depends on the grid level 
k, see below. We will use an inverse inequality to show that || • ||jv < CHI • |||fc holds. 
The convergence of the multigrid method (smoothing property and approximation 
property) will be shown in the norm ||| • |||fc. 

The main focus of the paper is the proof of the approximation property. Through¬ 
out the paper we will comment on the differences to the approach presented in |15| . 
which allows us to drop the requirement of solving (approximately) Poisson problems. 
The key idea of the proof of the approximation property follows classical approaches. 
However, the way of constructing the norms is non-standard. For simplicity, we follow 
the abstract framework introduced in m- First we introduce the following convenient 
notation. 

Notation 2.1. Throughout this paper, C > Q is a generic constant, independent 
of the grid level k and the choice of the parameter ft. For any scalars a and b, we 
write a<b (or b> a) if there is a constant C > 0 such that a < C b. We write b 
if a <b < a. 

Let the Hilbert spaces X, U and P (introduced above) be equipped with the 
following norms: 

II^IIa := IIKp)IIa := Ml + Ml, 

\W\\u •= ll'^llRqn) + /3|klli2(n) and 


Lemma 2.1 in m states the following stability result. 
Lemma 2.2. The relation 


II^IU< 


sup 

O^x^X 


B{x, x) 


<ll^llx 


(Al) 


holds for all X € X. 

Using the following notation, we can express the norms in a nicer way. 
Notation 2.3. For any Flilbert space A, the term A* denotes its dual space 
equipped with the dual norm 


\\u\\a* ■■= sup 




o^weA I|it||a 


where (u, •) := it(-) denotes the duality pairing. 

For any Flilbert space A and any scalar a > 0, the term a A denotes the space on 
the underlying set of the Hilbert space A equipped with the norm 

MlA--=a\\nrA. 


For any two Hilbert spaces A and B, the term Af] B denotes the space on the inter¬ 
section of the underlying sets, {u € AD B}, equipped with the norm 

MlnB ■■= IkllA + lklll 


and the term A B denotes the space on the algebraic sum of the underlying sets, 
{ui U 2 '■ ui G A, U 2 € B}, equipped with the norm 
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The spaces A*, a A, A D B and A + B are Hilbert spaces. The fact that A* 
is a Hilbert space follows directly from the Riesz representation theorem, see, e.g., 
Theorem 1.2 in [T]. The fact that a A is a Hilbert space is obvious and for the latter 
two see, e.g., Lemma 2.3.1 in 

We immediately see that the norm on U can be rewritten as follows: 

Ihllc/ = \\u\\H^{n)ni3L^{n)- 

To reformulate the norm on P, we need the following regularity assumption. 

(R) Regularity of the Stokes problem. Let / G and g G HQ{n) H Lo(^) 

arbitrary but fixed and {u,p) G [iLQ(fl)]^ x be the solution of the Stokes 

problem, i.e., such that 

(Vu, Vu)l2(o) + (p,Vu)L2(n) = (/,u)L2(f2) 

(Vu,p)L2(f2) = (5 ,p)l2(q) 

holds for all {u,p) G [7LQ(n)]^ x Lg(n). Then {u,p) G x H^{n) and 

Il'“llff2(n) + Ibliffi(f2) ^ Il/lli2(n) + ll5llHi(n)- 

Lemma 2.4. The regularity assumption (R) holds for 12 C being a convex 
polygonal domain. 

Proof. Theorem 2 in [10] states (in the notation of the present paper) that pro¬ 
vided / G [L^(12)]^, g G iL^(12) n Lo(^) ^ L^(f2) that 

lbllff2(o) -I- ||Vp|||2(n) ^ Il/lli2(n) + l|V5|||2(n) + lb ^5lli2(n); 

is satisfied, where (5 : 12 —>■ K is the distance to the closest corner of the polygonal 
domain 12. Lemma 2 in [10] states that |b~^ 5 ||L 2 (Q) < ||(?||iji(n) is satisfied for all 
g G 2 / 0 ( 12 ). Combining these results, we obtain that 

lblllr2(Q) -I- l!Vp|||2(n) < Il/Ili2(n) + llffll^i(n) 

is satisfied for all / G [L^(12)]^ and g G 2/o(12) H Lq(12). As p G Lq{D,) was assumed, 
Poincare’s inequality states that |b|l//i(n) ||Vp|b 2 (n), which finishes the proof. □ 

Note that we assume that g satisfies homogeneous Dirichlet boundary conditions. 
This condition can be weakened but it is not possible to drop such a condition com¬ 
pletely, cf. [To] . 

Lemma 2.5. On domains 12, where (R) is satisfied, 


lbl|p ~ \\p\\l^(Q)+0-'^H^(CI) 


holds for all p G Lq(12). 

For a proof of this lemma, see Theorem 3.2 in [1^ (this proof only needs (R) for 
/ e [L^(12)]^ and 5 = 0). 


The discretization of problem (2.11 is done using standard finite element tech¬ 
niques. We assume to have for fc = 0,1, 2,... a sequence of grids obtained by uniform 
refinement. The grid size on grid level k (length of the largest edge) is denoted by 
hfc. On each grid level k, we discretize the problem using the Galerkin approach, i.e., 
we have finite dimensional spaces Xk C X and consider the following problem. Find 
Xk G Xj- such that 


B{xk,Xk) = B{xk) lov all XkGXk. 


(2.3) 
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Following the concept of mixed finite elements, we choose Xk := Uk x Pk, where 
Uk C U and Pk C P are obtained by a proper discretization, see below. Using a nodal 
basis, we can represent this problem in matrix-vector notation as follows: 


Axfc =/^, 

where 


(2.4) 


Ak = 


Ak 

Bk 





and the matrices Ak and Bk represent the scalar products o(m, u) := (Vu, Xu)L^(n) + 
I3 {u,u)l2(q) and b{u,p) := (V • u,p)L 2 (n), respectively. 

Here and in what follows, any underlined quantity, like Xk, is the representation 
of the corresponding non-underlined quantity, here Xk, with respect to a nodal basis 
of the corresponding Hilbert space, here Xk- 

The next step is to show the discrete stability condition, i.e., that 


\\xk\\x < sup 


B{xk,Xk) 


^ W^kWx 


(Ala) 


holds for all Xk G Xk- 

To guarantee the discrete stability condition, we have to choose a discretization 
which is stable for the standard Stokes problem, particularly a discretization that 
satisfies the weak inf-sup condition, i.e.. 


(V • Uk,Pk)L'^{n) 

-n—n- 

o¥=ukeUk \\Uk\\L2{Q) 


^ l|Vpfc||L2(n) 


(2.5) 


should hold for all Pk & Pk- Note that this is a standard condition which guarantees 
that the chosen discretization is stable for the Stokes problem. In |5J it was shown 
that condition (2.5) is satisfied for the Taylor-Hood element (PI — P2-element) for 
polygonal domains where at least one vertex of each element is located in the interior 
of the domain. Here and in what follows we assume that the problem is discretized 
with the Taylor-Hood element and that the mesh satisfies the named condition. Using 
the weak-inf sup condition (2.5) we can show the following lemma. 

Lemma 2.6. If the problem is discretized using the Taylor-Hood element, condi¬ 
tion (Ala) is satisfied. 

Proof. The estimate (2.16) in [15] states that 

B{xk,Xk) 


\Xk\\Xk,k 


< 


sup 


holds for all Xk G Xk, where 

:= \\{uk,PkWxk,k ■■= 

■= ll^fcll^qn) +/3|l“fcllL2(n) and 

{Pk:^ ■ Wk)l2(^a) 

\\Pk\\p^k-= sup - ^72 -—- ^72 - . 

Here, the underlying sets of the Hilbert spaces Xi^k, Ui^k and Pi^k coincide with those 
of Xk, Uk and Pk, respectively. Observe || • ||( 7 i = || • |j( 7 . Lemma 2.2 in [15] states 

II ■ llfi fc ~ II ■ Ilf - This shows (Ala) as boundedness follows directly from (Al). □ 
Note that the grids are obtained by uniform refinement. So, the discrete subsets 
are nested, i.e., Uk C Uk-i-i, Pk U Pk+i and, consequently, Xk C Xk+i. 
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3. A coupled multigrid method. The problem (2.41 shall be solved using a 
multigrid method. Starting from an initial approximation one iterate of the 
multigrid method is given by the following two steps: 

• Smoothing procedure: Compute 


(0,m) 


(0,m-l) 


A\ 


" (4 


-Ak 




for TO = 1,..., z/ (3.1) 


with := The choice of the smoother (or, in other words, of the 

preconditioning matrix will be discussed below. 

• Coarse-grid correction: 

— Compute the defect f ^ — Ak 3±k''''^ restrict it to grid level k—\ using 


an restriction matrix I\ 


k-l. 


.^(1) ._ jk-l 

t-k-1 ■— -'fc 


— Determine the update by solving the coarse-grid problem 

Ak-i ^k-i — ti^k-i (3-2) 


approximately. 

— Prolongate to the grid level k using an prolongation matrix I^_i 
and add the result to the previous iterate: 


„(i) 


„(o.C 


^k 


^k 


rk (1) 

■ ^k-1 S.k-1- 


As we have assumed to have nested spaces, the intergrid-transfer matrices Ik-i and 
are chosen in a canonical way: Ik-i is the canonical embedding and the restric¬ 
tion is its transposed, i.e., = (/|_i)^. 

If the problem on the coarser grid is solved exactly (two-grid method), the coarse- 
grid correction is given by 


.(1) 




-I- A~'^ T‘ 

\ -1 I AAu 1 -L ) 


-1 




-A. 


soy, 


k Afc 


(3.3) 


In practice, the problem ( |3.2[ ) is approximately solved by applying one step (V-cycle) 
or two steps (W-cycle) of the multigrid method, recursively. On grid level fc = 0, the 
problem (3.21 is solved exactly. 


To construct a multigrid convergence result based on Hackbusch’s splitting of the 
analysis into smoothing property and approximation property, we have to introduce 
an appropriate framework. 

Convergence is shown in the following L^-like norms ||| • |||fc: 


lll^^fclllfc ■ ll3iA:ll£)j ■ i^kXki !S-k)i^ 1 


where 


:= 


(hk ^ + P)Mu,k 


(/^ + ) ^^p,k 


(3.4) 


the matrices Mjj^k and Mp^k are mass-matrices, representing the L^-inner product in 
Uk and Pk, respectively, and (•, •)^2 is the Euclidean scalar product. 

The smoothing property and the approximation property, which we will show 
below, read as follows. 
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Smoothing property: 


B 


sup ■ 


(4°'“ 


- X^,Xk 


IPfclllA 


holds for some function t](i') with limj^ 


< ?7(f^)lll4°^-4lllfc (3-5) 

rj{v) = 0. Here and in what follows, 


x*). € Xk is the exact solution of the discretized problem (2.3). 
• Approximation property: 


B 


( 4 °’ 




- Xl,Xk 


IFfcllfc 


(3.6) 


\\\x^k^ - xl\\\k < Ca sup ■ 

Xk^Xk 

holds for some constant Ca > 0. 

If we combine both conditions, we obtain 

Ill4^^ - 4lllfe < qH\\\Xk^ - xl\\\k, 

where q{v) = CAVi^)- ^ large enough, we obtain q{i’) < 1, i.e., the convergence 
of the two-grid method. The convergence of the W-cycle multigrid method can be 
shown under mild assumptions, see e.g. [3]. 

Remark 3.1. The norm ||| • |||fe is the discrete analog of X^^k ■= x Po^k ■= 
Uk X Pk, equipped with norms 


Ixkfxo.k ■■= IIK,Pfc)ll4.. := 11^^114. + Ib^llPc,.’ 

|2 


WWuo, 


■— lbfellL2(o) +/?||Mfe||L2(Q) and 
{PkiWkYkif^Q^ 


\\Pk\\pok'= sup 

0^Wk&L^{fi) "'k 


(3.7) 


hk^WwkP 


L2(n) 

k w^) ■^Ibfc||i2(n). 


+ /3 Wk 




= h-,\h:^ ' 


Note that this norm is obtained from the norm || • ||x by “replacing all differentials by 
^-ij; jg g common construction principle, cf. m- However, in m the norm 


\xk\\l^,,--= \\iuk,Pk)\\l,^ 


WPkWPiA ■= sup 


:= Wukfu^^k + lMp., 
b ^Ibfclli2(n) 

(Pfe,V • Wfc)i2(o) 


(3.8) 


o^wk&Uk ll'*Pfcll/ri(n) + Pll'^fellL2('Q) 

was used. Note that the norm for the velocity field coincides with our choice but the 
norm for the pressure distribution is still the original norm, as introduced in (2.6). 

4. Smoother and proof of the smoothing property. The choice of an ap¬ 
propriate smoother is a key issue in constructing a coupled multigrid method for an 
indefinite problem. In this paper, we introduce two kinds of smoothers. The first 
smoother is appropriate for a large class of problems including the model problem: 
the normal equation smoother, cf. |7], which reads as follows. 


(0,1 

^k 


_ 

Tk 


+ T Ufc/ifc ^ (/^ - Ak 4 °’’" 




for TO = 1,... ,p. 
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Here, a fixed r > 0 has to be chosen such that the spectral radius p{TA^^Ak) is 
bounded away from 2 on all grid levels k and for all choices of the parameter /3. 

Remark 4.1. In Section 6.1 in m a normal equation smoother has been pro¬ 
posed. (There, this kind of smoothers was called distributive smoother, j Note that 
this class of smoothers depends on the choice of the Hilbert space norm. In the present 
paper, a scaled if-norm is used. Therefore, the Riesz isomorphism Ck can be easily 
inverted. However, in m a different norm was chosen, cf. Remark \3.1\ For that 


norm, the realization of the inverse of the Riesz isomorphism involves the solution of 
one Poisson problem (for the pressure variable). So, two such problems have to be 
solved (approximately) in each smoothing step, cf. Section 8 in ]15^ . 

Using a standard inverse inequality, one can show that 


\Xk\\x 


< 


(A 2 ) 


is satisfied for all Xk € Xk. Based on this result, using an eigenvalue analysis one can 
show the following lemma, cf. |7]. 

Lemma 4.2. The damping parameter r > 0 can be chosen independently of the 
grid level k and the choice of the parameter (3 such that 

Tp{Af^Ak) < 2 - e < 2, 

holds for some constant e > 0. For this choice of t, there is a constant Cs > 0, inde¬ 
pendent of the grid level k and the choice of the parameter f, such that the smoothing 
property (3.5) is satisfied with rate rjiy) := Csv~^l'^. 


Certainly, the smoothing procedure (3.5) should be efficient-to-apply. Using the 


fact, that the mass matrices Mu^k and Mp^k in (3.4) and their diagonals are spec¬ 
trally equivalent under weak assumptions, for the practical realization of the normal 
equation smoother these mass matrices can be replaced by their diagonals. 

The second smoother, which we propose, is a Uzawa type smoother, cf. m- Here, 
one step of the smoother to compute = iu^k ’^\pk 

reads as follows: 


(0,m) _ ,JO,m) = 


(0,m-l/2) 




:= af + rAp (4 - 


- '"4-' (a. - 


(0,m)_(0,m—1) 




■ rife 


+ ^ (/fe “ ^kU'k’ 


m—1) 


T (0,m) 


^kPk 


). 


where Ak and Sk are the (l,l)-block and the (2,2)-block of Ck, respectively. The 
damping parameters r and a have to be chosen independently of the choice of fi 


compact notation (3.1), where 


and of the grid level k such that (4.1) holds. The smoother can be rewritten in the 


Ak := 


T ^Ak 
Bk 


TBkAr^B^ - a-^Sk 


The smoothing property can be shown using the theory introduced in m- 

Lemma 4.3. Let Ak and Sk be the (l,l)-block and the (2,2)-block of Ck, respec¬ 
tively. Then r > 0 and cr > 0 can be chosen independently of the grid level and the 
choice of the parameter fi such that 


-1 


Ak A Ak 


and 


a Sk > rBkA^ 


Bl- 


(4.1) 
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For this choice of t and a, there is a constant Cs > 0, independent of the grid 
level k and the choice of the parameter (3, such that the smoothing property (3.5) is 
satisfied with rate r]{v) '■= 


Proof. The fact that r > 0 and a > 0 can be chosen independently of the grid 
level and the choice of the parameter /3 follows from standard inverse inequalities. 
For the smoothing property, we apply Theorem 4 in m with the choice JCk '■= 
and /Sfe := ■ This immediately implies 


As = B{xk,Xk) < \\xk\\x\\xk\\x < IllxfcllUpfclllfc = \\xj,\\cAx^\\ck holds for 

all Xk,Xk € Xk, we obtain AkCff^^'^\e‘ < 1, which finishes the proof. □ 

Remark 4.4. Uzawa type smoothers have also been considered in Section 6.2 
in m- As for the case of normal equation smoothers (distributive smoothers) due to 
the particular norm that was chosen in m, again a discrete Poisson problem for the 
pressure variable has to be solved in each smoothing step, cf. Section 8 in m- 

5. Proof of the approximation property. The proof of the approximation 
property ( |3.6[ ) is done using the approximation theorem introduced in m which 
requires besides the conditions (Al) and (Ala) two more conditions (conditions (A3) 
and (A4)) involving, besides the Hilbert space X, two more Hilbert spaces X-^k ■= 
{X_, II • ||jf_ A and X^^k ■= {X+, || • ||x+ which are chosen as follows. 

As weaker space we choose A_ := [/_ x P_, where U- := [£^(12)]^ and P_ := 
[71g(r2) n £ 0 ( 12 )]*. These Hilbert spaces are equipped with norms 

PIIL. := II(«,£)IIL.. := Pllu_.. + lbllp-..^ 

Ibllu_,fc ■= and 

■= ^fe^lbllpi-i(n)]*-i-/3-iL2(n)- 


Remark 5.1. The idea behind the construetion of the norm || • ||j(:_ ^ is is to take 
the norm || • ||x vLnd “replace all occurrences of by h^^L^ and all occurrences of 
by h~j(^H~^”. This is different to the idea of constructing the norm ||| • |||fe, where 
“only the -terms are replaced by h^^L^ and the -terms are kept unchanged”. 

In the following, we show the approximation property in the norm || • ||x_ i.e.. 


£ 1 ) 






-Xfcllx.,*, < C'a sup . 


-Xk,Xk, 


\Xk\\x.., 


(5.1) 


We will show below that this version of the approximation property is stronger than 
the required estimate (3.6), cf. Lemma 5.9 In classical approaches, cf. |20l I16| . the 


approximation property is directly shown in the norm ||| • |||fe. 

Note that dual spaces are (A_)* := (U-)* x {P-)*t where (U-)* = [A^(f2)]^ and 
(P_)* = iJg(r2) n Lq(J1), equipped with norms 


ll/ll((7_.fe)* = ^fcll/llLRa)+/3-ipi(G) 
lbll(p_,fe)* = ^fclbliffi(n)n/3Lg(n)' 
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As stronger space we choose A+ := U+ x P+, where C/+ := [i?^(f2)ni?g(r2)]^ and 
P+--= n Pg(ll), equipped with norms 

ll^llx++ := ll(«,P)llx++ := hll&++ + lbllp++, 

lbllp++ ■= and 

lbllp++ ■= ^fclblipi(n)+/3-ip2(o)- 

For showing the approximation property, we need some auxiliary results. The first 
result is a standard approximation error estimate: 

Theorem 5.2. On all grid levels k, the approximation error estimate 

inf lb - Xfellx < |b||x+fc for all X G X+ (A3) 

XfcGAfc 


is satisfied. 

The proof of this theorem is rather easy, as it can be shown for U and P separately. 
For completeness, we give a proof of this theorem in the Appendix. The next step is 
to show a regularity result for the generalized Stokes problem: 

(A4) For all grid levels k and all T G (A_)*, the solution of the problem, 

find xjr G X such that B{xj^,x) = P{x) for all x G X, (5.2) 


satisfies xjr g X+ and |bjr||x++ < ll-^^ll(x_+)*• 

For showing (A4), we need a standard regularity result for the Poisson problem. 
(Rl) Regularity of the Poisson problem. For all g G the solution of the 

problem, find pg G n Lq{Q) such that 

(Vpg, Vp)^i(n) = {g,p)L^(n) for all p e H^{n)nLl{n), 

satisfies pg G and |bgl|p2(n) b hh^in)- 

Lemma 5.3. The regularity assumption (Rl) holds for 12 C being a convex 
polygonal domain. 

For a proof oft this lemma, see, e.g., [S]. The first step of the proof of (A4), is 
to show that xjr is in the desired space. 

Lemma 5.4. Suppose that 12 is such that the regularity assumption (R) holds. 
Let T G {X-)* and let xjr be the solution of ( |5.2[ ). Then xjr G A+. 

Proof Let P{u,p) := (/,u)l 2 (o) + {g,p)L^Q) for / G [L‘^{0.)]'^ and g G H^{n) n 
Lfj{Tl). Rewrite the problem as follows: 


(VujF, V£t)i2(n) + (pp, V • u)l2(q) = (/, ■u)l2(q) —/?(up, {t)i2(Q) 

(V ■ up,p)L2(n) = {g,p)L^(Q) 


for all u G U and p G P. As / — /3up G [L^(12)]^ and g G ilo(12) H Lg(12), we obtain 
using regularity assumption (R) immediately that xp G A+. □ 

Note that the combination of the argument used in the proof of this lemma 
and condition (Al) implies ||a:p|jp_^_j^ < C'(/3)|| where C{f3) is a constant 

depending on /3. For showing a robust estimate, we need to do some more work. 

Lemma 5.5. Suppose that 12 is such that the assumptions (R) and (Rl) are 
satisfied. Let T € and let xjr = (ujr.,‘pjr) he the solution of (5.2). Then 
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Lemma 5.6. Suppose that O is such that the assumption (R) is satisfied. Let 


J- S (^-)* and let xj^ = (ujr,pjr) be the solution of (5.2|. Then 


\WM\u+_k - II•^II(X_,0•I|2;J=■||X4. 


The proofs of these two lemmas can be found in the Appendix. Using them, we 
can show (A4): 

Theorem 5.7. If LI is such that the regularity assumptions (R) and (Rl) are 
satisfied, then the regularity statement (A4) holds. 

Proof. Note that Lemma [5.4| already states that xjr g A+. It remains to show 
||a;j='||A:+ k ^ 11-^11 (a:_ t,)* Note that Lemma 5.5 implies immediately 




= Mi;,,, + llPrWK.. ^ Mi;,,, + 

If we combine this result with the statement of Lemma 15.61 we obtain 




for some constant C > 0 (independent of k and /3) which implies 


<i(c+\/4C + C2) ||J-||p 


and the desired estimate. This finishes the proof. □ 

As we have shown that the conditions (Al), (Ala), (A3) and (A4) are satisfied, 
we can apply Theorem 4.1 in |20| and obtain: 

Theorem 5.8. Assume that LI is such that (R) and (Rl) are satisfied and as¬ 
sume that the problem is discretized using the Taylor-Hood element. Then the coarse- 
grid correction (3.31 satisfies the approximation property ([53 where Ca only depends 
on the constants that appear in the conditions (Al), (Ala), (A3) and (A4). 

This shows (5.11, the approximation property in the norm |j • however we 

need (3.6 1 , the approximation property in the norm ||| • |||fc. The final step is to relate 


these two norms: 

Lemma 5.9. The estimate 




< 


\Xk\\x_,k 


(5.3) 


holds for all Xk G Xk 

Proof. First note that it suffices to show 
||Pfc||p_ k- obtain using the definition of || • 
and a standard inverse inequality 


u,k < \\uk\\u-k and 


ip.fc 

|c/_ the choice u := Uk e [H^{Lt)]'^ 


\\Uk\\u. , = 


sup 


iuk,ufi 


L2(0) 


(wfe, Rfc)i2(f2) 

inf^GK(n)]= hl\\w - Uk\\l2(^a) + ^fcll^ll^-ipi(n) 


> 


{uk,Uk^'^ 


L2(n) 


hl\W'^ 


L2(a)+/3-ipi(n) 


> 


{uk, Uk, 


L2(n) 


inf. 


ujfcG 


2 


> 


{Uk,Uk)\2(Jl'^ 


infujfcGPfc '*^felli2(f2) + ll'a’fe|l^-iL2(Q) 


— (/5 + ■ 


Ffe|lL2(0) 


2 

U.k- 
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For the proof of the inequality for the pressure variable, we define on each grid 
level a function —>■ M as follows. The function -0^ is 0 on all vertices which 

are located on dfl, 1 on all other vertices (interior vertices) and linear within each 
element (Courant element). 

In what follows, Pk is the space of functions obtained by discretizing the space 
H^{n) using the standard Courant element. Note that ~ as we have used the Taylor 
Hood element - the identity Pk = {pk + a : Pk & Pk and a € K} holds. 

Note that we obtain using the definition of || • ||p_ and q C To(H) that 


\\Pk\\l_,= 


sup 


{pk,qf 


L2(a) 


sup 


(pfc -a,9)i2(n) 


o#96ff(i(n)nL2(n) ^felkllp-i(n)n,SL2(n) ^k\\q\\ 

is satisfied, where a G K is chosen such that {pk — a)ipk G Poi^)- Observe that 
q := [pk — a)'ijjk G iFo(H) C ^^(H). Using this choice, we obtain 


\\Pk\\ 


^ > 
P-,k - 


{Pk - a, {pk - a)'tjjk)l 2 (^Q-^ 


hlWiPk - a)0/c||^i(a) + ^hl\\{pk - a)iJk\\l2(^^^' 


Using \ipk\ < 1, we obtain \\{pk - a)0fe||^2(Q) < \\Pk - a|li 2 (Q). The following two 
inequalities hold: 

llV'fcPfcllpi(n) ^ ^lbfcllL 2 (n) and (5.4) 

iPk,'4’kPk)L^n) > {Pk,Pk)L^Q)- (5.5) 

Proofs are given in the Appendix. Using these inequalities and pk G Lg(U), we obtain 

(Pfc-a,Pfc-a)|2(n) 


\\Pk\\ 


^ ^ Ibfc - a|||2(n) + I3hl\\pk - aWhi^n) 

= {I + f3hl) ^\\pk — a\\'^2(^Q) > {I + (3hk) ^|bfc||i2(Q) = |||p/c|||p_fc. 


which finishes the proof. □ 


Using this lemma and the approximation property (5.11, we obtain the approxi¬ 


mation property (3.6). Together with the smoothing property, we obtain: 


Theorem 5.10. Assume that 

• n is such that the regularity assumptions (R) and (Rl) are satisfied, 

• the problem is discretized using the Taylor-Hood element and 

• one of the smoothers proposed in this paper is used. 

Then 


Il4^^ -4lllfc < qip) 


00) 


holds with q{v) := CsCav~^^'^, where the constants Ca and Cs are independent of 
the grid level k and the choice of the parameter [3. For v > (Cg Ca)"^, the two-grid 
method converges with convergence rate q{v) < 1. 

The convergence of the W-cycle multigrid follows under weak assumptions, cf. |1]. 


6. Numerical Results. In this section, we illustrate the convergence theory 
presented within this paper with numerical results. For the numerical experiments, the 
domain U was chosen to be the unit square fl := (0,1)^. As mentioned in Section]^ the 
weak inf-sup-condition (2.51 can be shown for the Taylor-Hood element only if at least 
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one vertex of each element is in the interior of the domain O. As this is not satisfied 
for the the standard decomposition of the unit square into two triangular elements, 
we choose the coarsest grid level fc = 0 to be a decomposition of the domain O into 
8 triangles, as seen in Figure |6.1[ The grid levels k = 1,2,... were constructed by 
uniform refinement, i.e., every triangle was decomposed into four subtriangles. 



Figure 6.1. Discretization on grid levels k = 1 and k = 2, where the squares denote the degrees 
of freedom of (the components of) the velocity field u and the dots denote the degrees of freedom of 
pressure distribution p 


The right-hand-side functions / and g have been chosen such that the solution of 
the problem on grid level k is the L^-projection of the exact solution 

■= 4>{Ci,^2) Cl ) ^(^1)^2) := ^(>(^1)6), 

into Afc, where (/)(Ci,C 2 ) := max |o, min |l, 2 - + (C 2 - 

For solving the discretized problem, we have used the proposed W-cycle multigrid 
method. As mentioned above, the matrix Ck, introduced in Section|^can be replaced 
by any spectrally equivalent matrix. As Ck should be easy-to-invert, Ck should be a 
diagonal matrix. The convergence theory indicates how Ck should be scaled, however 
only up to constants. Numerical experiments indicated that it is reasonable to choose 

^ ^ ^ , where Ak := diag Au and Sk ■= diag BkA'^'^Bl, (6.1) 


i.e., to choose these constants as they appear naturally in the diagonal of Ak, the (1,1)- 
block of Ak, and in the diagonal of the inexact Schur complement. As this choice of 
Ck is spectrally equivalent to the matrix Ck, introduced in Section this choice is 
still covered by the convergence theory. For this choice, the damping parameter was 
chosen to be r = 0.35, for all grid levels k and all choices of (3. 

For the Uzawa type smoother, the matrices Ak and Sk have been chosen as 


introduced in (6.1) and the damping parameters have been chosen to be r = cr = 0.8 


for all grid levels k and all choices of j3. Again, the choice of this section is covered 
by the convergence theory. 

The number of iterations and the convergence rate were measured as follows: we 
start with and measure the reduction of the error in each step using the norm 
III • lllfe. The iteration was stopped when the initial error was reduced by a factor of 
e = 10“®. The convergence rate q is the mean convergence rate in this iteration, i.e.. 




l/n 


<? = 




.( 0 ) 


- XI 


where n is the number of iterations needed to reach the stopping criterion. Here, 
is the exact solution and x). is the i-th iterate. 
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V = 1 + \ V = 2 2 

u = 

3-f 3 

1/ = 4 -1- 4 

z/ = 8 -k 8 

ly = 

16-k 16 

n q n q 

n 

q 

n q 

n q 

n 

q 

Normal equation smoother 

88 0.789 46 0.634 

30 

0.496 

24 0.412 

17 0.293 

11 

0.148 

Uzawa type smoother 

divergent divergent 

14 

0.212 

20 0.347 

6 0.031 

5 

0.008 

Asymptotic decay of the convergence rate predicted by theory (z/ ^ 

/2) 


0.707 0.500 


0.408 

0.353 

0.250 


0.177 


Table 6.1 

Number of iterations n and mean convergence rate q for the normal equation smoother and the 
Uzawa type smoother depending with v = v-pre + ^'post smoothing steps on grid level k = 4 for /3 = 1 




/3 = 

n 

0 

q 

/3 = 

n 

10^ 

q 

/3 = 

n 

10^ 

q 

/3 = 

n 

10^ 

q 

P = 

n 

10® 

q 

f3 = 

n 

IQio 

q 

k 

= 4 

30 

0.496 

30 

0.493 

30 

0.496 

68 

0.736 

71 

0.745 

71 

0.745 

k 

= 5 

29 

0.489 

29 

0.488 

22 

0.388 

64 

0.722 

70 

0.744 

71 

0.745 

k 

= 6 

29 

0.484 

29 

0.486 

27 

0.457 

52 

0.670 

70 

0.743 

71 

0.745 

k 

= 7 

28 

0.475 

28 

0.475 

28 

0.470 

35 

0.553 

70 

0.742 

71 

0.745 

k 

= 8 

28 

0.469 

28 

0.469 

28 

0.468 

20 

0.347 

67 

0.732 

71 

0.746 


Table 6.2 

Number of iterations n and mean convergence rate q for the normal equation smoother with 
= 3 + 3 smoothing steps 


In Table 6.1 we compare for a fixed grid level (level fc = 4) and the fixed choice 
/3 = 1 the convergence rates for several choices of the number of pre- and post¬ 
smoothing steps. We see that the convergence rate behaves approximately like 
if the number of smoothing steps is increased. This is consistent with the theory which 
guarantees the convergence rate being bounded by We observe that the 

preconditioned normal equation smoother already converges for v = 1 1 smoothing 

steps, while for the Uzawa type smoother z/ = 3 -I- 3 smoothing steps are necessary. 

In Tables 6.2 and 6.3 we compare various grid levels k and choices of the param¬ 
eter fi. We have used a fixed choice of z/ = 3 -I- 3 smoothing steps. First we observe 
that, for both smoothers, the number of iterations seems to be well-bounded for all 
grid levels k which yields an optimal convergence behavior. Moreover, we see that the 
number of iterations is also well-bounded for a wide range of choices of the parame¬ 
ter /3, i.e., we observe also robust convergence as predicted by the convergence theory. 
Comparing both kinds of smoothers, we see that the Uzawa type smoother leads to 
much faster convergence rates than the preconditioned normal equation smoother. 
Note that moreover the computational complexity of the Uzawa type smoother (per 
iteration) is slightly smaller than the complexity of the normal equation smoother. 

It has to be mentioned that for the model problem, also the V-cycle multigrid 
method converges with rates comparable to the convergence rates of the W-cycle 
multigrid method. However, the V-cycle is not covered by the convergence theory. 
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/? = 
n 

0 

q 

n 

10^ 

q 

13 = 

n 

10-* 

q 

/3 

n 

= 10® 

q 

13 

n 

= 10® 

q 

P 

n 

= 1010 

q 

k 

= 4 

14 

0.212 

13 

0.201 

6 

0.030 

7 

0.051 

8 

0.059 

8 

0.059 

k 

= 5 

13 

0.194 

13 

0.193 

9 

0.095 

7 

0.043 

7 

0.050 

7 

0.050 

k 

= 6 

12 

0.176 

12 

0.176 

11 

0.145 

6 

0.024 

7 

0.047 

7 

0.048 

k 

= 7 

12 

0.166 

12 

0.166 

11 

0.150 

5 

0.013 

7 

0.044 

7 

0.045 

k 

= 8 

11 

0.147 

11 

0.147 

11 

0.145 

8 

0.058 

7 

0.038 

7 

0.043 


Table 6.3 

Number of iterations n and mean convergence rate q for the Uzawa type smoother with i/ = 3 + 3 
smoothing steps 


The numerical experiments done by the author have shown that the convergence 
rates can be improved slightly by adjusting the choice of the parameters to the grid 
levels and the choice of /3. However, the main goal of this paper is to show that the 
proposed method also works well for fixed choices of the parameter. 


7. Conclusions and Further Work. In the present paper we have proposed 
a coupled multigrid solver for the generalized Stokes problem where the smoothing 
property is shown in the scaled L^-norm ||| • |||fe. This allows to construct a multigrid 
method where the smoother is a simple linear iteration (which consists only of divisions 
and the multiplication of vectors with the system matrix Ak)- In the present paper, 
a preconditioned normal equation smoother and an Uzawa type smoother have been 
chosen but it seems possible to find also other smoothers which satisfy the smoothing 
property in the norm ||| • |||fc. 

The convergence rates observed for the multigrid method proposed in the present 
paper are comparable with the rates observed for the methods proposed in [15] . Note 
that for applying the methods proposed in the named paper, it is necessary to solve a 
Poisson problem in each smoothing step. This is not needed for the method proposed 
in the present paper. The main contribution of this paper is a new way of setting up 
the norms where the convergence is sown in. The technique of the convergence proof 
that has been applied in the present paper is also extendable to the Stokes control 
problem, cf. m- 

Acknowledgments. The author thanks Markus Kollmann for providing code 
used to compute the numerical results presented in this paper. Moreover, the support 
of the Austrian Science Funds (FWF), the University of Oxford and the Technische 
Universitat Chemnitz are gratefully acknowledged. The author thanks the referee for 
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8. Appendix. Proof of Theorem 5.2. Note that it suffices to show approxima¬ 


tion error results for the individual variables separately. Using a standard interpola¬ 
tion operator H^ : —>■ U^, we obtain 


||u — nfcu|||2(Q) < ^fcl|w||ffi(n) and ||m — Hfeulj 


for all u € [i7^(r2)]^ and therefore 

inf ||u - UkWl < l|w - ^kuWl = ||m - + l^\W - nfcM|li 2 (n) 
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A similar argument is possible for the pressure distribution. The standard approxi¬ 
mation error estimates 


plek ^ and ^in^^ Ib-Pfcllffi(n) < ^fclbllff 2 (o) 

hold for all p G H^(Q) or p G respectively. They imply 

inf |b-Pfe||p= inf inf Ibi-Pi,fc|li 2 (n)+ inf |b 2 - P 2 .fc|||-ipi(n) 

PkGPk pieL^(n) pi.kSPk ^ ’ P2,kePk ^ ^ ' 

P2eH\n) 

Pl+P2=P 

< inf inf Ibi-Pi.fc|li 2 (n)+ inf Ib 2 -P 2 .fe||^-ipi(n) 

pieH^in)pi,kePk '• ' P2,kePk ^ ^ ' 

P2GH^(n) 

Pl+P2=P 

^ inf Ibillpi(n) + Ib2||^-ip2(n) = fifclbllHi(n)+/3-iff2(n), 

pi&H (n) 

P2GH^{n) 

Pl+P2=P 


for all p G n Ln(n), which finishes the proof. □ 


5.5 


Proof of Lemma 
g G n LQ(r2).' Clioose /2 

[L^(r2)]^. Observe that 


Let J^{u,p) := {f,u)L^(n) + {g,p)mn) for / G [T^(0)]2 and 
G [iLg(O)]^ arbitrarily and define fi := f — f2 G 


{pjr, V • u)l2(CI) = -(Vup, Vu)p 2 (Q) - /3{ujr, u)p 2 (Q) -|- (/l -|- f2 , u) 1,2 

is satisfied for all u G [iLQ(O)]^. This is equivalent to 

-(Vpp, u)l 2 (q) = (Aujr, u)p, 2 (^q) - u)l 2 (q) -I- (/l -1- / 2 , u)l^(Q,) 

for all u G [iLg(O)]^. Note that Vpp, Aup, up, /i and /2 are in [Lf(Vl)Y. Therefore, 
the above statement holds for all u G [L^(0)]^ (because [iLg (0)]^ is dense in [L^(0)]^), 
particularly for u := Vp, where p G Lf^(fl) n Lq(0). So, we obtain 

- (VPP, Vp)L2(n) = (Attp, Vp)p2(Q) - fl(ujr, Vp)p2(Q) -I- (/l -I- /2, Vp)p2(Q) (8.1) 

for all p G JL^(fi) fl Lg(r2). 

Let Pi G H^(Ll) n Lq( 0) be such that 

- (Vpi, Vp)p2(f2) = (Alip, Vp)p2(Q) -I- (/i, Vp)p2(f2) (8.2) 

holds for all p G H^{Ll) n Lg(0). 

Note that Auj^ + fi G [L^(r2)]^ and therefore the right-hand-side is a functional in 
[7Lg(0)]*. So, existence and uniqueness of pi G H^{Lt) n Lg(0) is guaranteed. Using 
the choice p := pi, we obtain 


l|Vpi|b2(n) < ||AMp||i2(o) -I- ||/i||l2(o) < \\ujr\\H2(n) + ||/i||L2(n)- 
Using Poincare’s inequality, we obtain further 

IbillffTfi) ^ lb.7=-|lr/2(Q) -I- ||/l|b2(Q). 

Let p2 G H^{Ll) n Lg(U) be such that 

- (Vp2, Vp)p2(Q) = -/?(up, Vp)i2(n) + (/2, Vp)i2(Q) 


(8.3) 
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for all p G 77^(0) n This implies, as ujr g [77g(0)]^ and /2 g [77 q( 0)]^, that 


- (Vp2, Vp)l2(q) = /3(V • MjF,p)L2(a) - (V • /2 ,p)l2(o) (8.4) 

holds. As /3V • ujr — W ■ f 2 £ L‘^{Vl), existence and uniqueness of p 2 is guaranteed. 
Condition (Rl) implies moreover p 2 G and 


lb2||ff2(Q) < /3^||V • u^||i2(Q) + ||V • /2||i2(n) < + ll/all^qn)- 


Note that from (8.2) and (^8.4|, we obtain 


(8.5) 


(V(pi +P 2 ),Vp)l 2 (q) = (Vpjr,Vp)l^ 2 (Q) 


is satisfied for all p G 77^(0)nLg(r2), which implies (because pj^ G 77q( 0) and pi +p 2 G 
Lg(0)) that pjr =p^-\-p 2 is satisfied. So, we have using (8.3) and (|8.5[) 


< p 4?+92 ^fell9illRi(£2)+^fcll92|||-i//2(o) 

9iGRbS^)nL’o(fI) 

92GR^(a)nL^(n) 

^ ^fclb.Tll/f-2(f2) + /lfe|bjG4ffi(0) + ^fcll/l|lL2(Q) + /lfe||/24-iRi(S2) 

= ^fclb.T|lb2(n)n/3Ri(n) + ^fcll/il!i2(n) + ^fell/2|||-iRbo)' 


As /2 £ [77q( 0)]^ was chosen arbitrarily, we can take the infimum over all / 2 , which 
finishes the proof. □ 

Proof of Lemma 5.6 As 770(0) is dense in L^(0), for u G [77^(0)]^ the func¬ 
tion —Au G [L^{Ll)Y can be approximated by some function w'^ G [77q( 0)]^ such 
that ||(—A)u — w'^||| 2 (Q) < e. So, we can introduce an operator —A'^ : [77^(0)]^ —)■ 
[77g(0)]^, which approximates — Au in [77 q( 0)]^, satisfying the error estimate 


(-A)u- (-A")u||i 2 (o) < 


< e 


for all u G So, —A^u approaches — Au in the L^-sense for e — > 0. (As 

—Au ^ [77g(0)]^ in general, |j — A'^u||//i(q) is not convergent for e —t 0.) 

Analogously, we introduce the operator : 77^(0) —>■ [77g(0)]^, which approxi¬ 
mates Vp in [77g(0)]^, satisfying the error estimate 

l|Vp-V'p||^2(n) < e 


for all p G 77^ (n). So, V'^p approaches Vp in the L^-sense for e —?> 0. (As Vp ^ 
[77g(n)]^ in general, || V'^p||//i(n) is not convergent for e —>■ 0.) 

The idea of this proof is to show that for all e > 0 there is some x‘^ G X such that 

xY > - e(l + and (8.6) 

< hfmix_,,]4xx\\x^,, + e(l + P^^4h4m\(x_,,y + e' (8.7) 


holds. As B{xjr,xY = F{x4, the statement of the Lemma follows for e —)■ 0. 
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In the following, we show (8.6) and (8.7) for ;= (—A'ujr, V • First, we 

show (8.6). Here, we estimate the summands of B{xj^,x'^) separately. We have 


(Vujr, = (Aujr, A'^Mjr)j;,2(f2) 

> (Aujr, Aujr)L 2 (n) - e|| Aujr||L 2 (n) > ||ujr||^ 2 (o) - 


where we use for the first step that A*^ maps into Hq{VI). For the second step, we use 
the Cauchy-Schwarz inequality and the upper bound for || — Aujr — (—A'^)uj 2 ||L 2 (f 2 p 
Friedrichs’ inequality is used for the third step. Using similar arguments, we obtain 
for the next summand 


> fiiujr, —AujF)L2(n) — e/311^x11^2(0) 

= - e/3||ujr||L2(n) > /3||wjF||^i(n) “ e/3/i^^||a::jr||x+,fc- 

We obtain for the last two summands, corresponding to the (2,1)- and the (l,2)-block, 
again using similar arguments, including the bound for ||Vpx ~ V'^px||L 2 (n), that 

(V • ux, V • V'^px)L2(n) -I- (V • (—A'^)ux,px)L2(n) 

= -(VV • ux, V'^px)L2(n) + (A'^ux, Vpx)L2(n) 

> -(VV • ujr, V'^px)L2(n) + (A'^ux, V'^px)L 2 (n) - e|| A^ux||L 2 (f 2 ) 

holds. Using integration by parts, the Cauchy-Schwarz inequality and the upper 
bound for || — Aux — (~A'^)ux||L 2 (n), we obtain further 

(V • uxj V • V'^px)L2(a) -I- (V • (-A'^)mxiPj2)j;,2(q) 

> (Vux, VV'^px)l 2 -I- (Aux, V'^px)l 2 -e(||AMx||L 2 -I- ||V'^px||l 2 -I-e) 

= (Vux, VV'^px)l 2 - (Vux, VV’^px)l 2 -e(||Aux||L 2 + ||V‘^px||l 2 -l-e). 

= 0 

Finally, we obtain using the upper bound for ||Vpx ~ e^nd the definition 

of II • IIx+,^ that 


(V • MX, V • V'px)L2(a) -I- (V • (-A'^)Mx,px)L2(n) 

> -e(||AMx||L2(n) + ||Vpx||L2(n) + e) > -e(||Mx||x2(Q) -|- Ibxll^qn) + e) 

> -e(l -f /3^^^)h^^||a:x||A+,fc - e^ 


holds. 

(/, M)L2(n) 


So, we have shown (|8.6|. The next step is to show (8.7). Let J-{u,p) := 
X2(U)]2 and g G H^i^) C L§(U). Let /a € [H^in)]^ 


(5,.p)L2(n) for / e 
and fi'=f — / 2 - Using the same arguments as above. 


(/l, -A'^Mx)l 2 (q) < ||/l||L2(Q)||ux||ff2(Q) -I- e||/l||L2(Q) 

holds as well as 


(/ 2 , — A'^Mx)l 2 (q) < (V/ 2 , VMx)L 2 (n) + e||/ 2 ||L 2 (f 2 ) 

^ ll/ 2 ||/ 3 -ir/i(n)l|wx||/ 3 jri(n) + e^^'^^ll/ 2 ||/ 3 - 

This implies 

(/, -A'^Mx)L2(n) 

^ ll/llL2(a)+/3-iffi(n)l|wxllH2(Q)n^xi(n) + e(l + /3^'^^)||/||L2(Q)+^-ij/i(n) 
^ ll/llL2(a)+/3-ifl-i(n)ll'Wx||H2(o)n/3i/i(n) + + P^^'^)\\^\\[x_^kY■ 
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Let p 2 G and pi ■— p — p 2 . We have 

(gi, V • V'^pi)L 2 (n) = -(V 5 , < ||g||_H-i(n)|bi|l/i-i(n) + e||5||Hi(n)- 

Using g G Hq(^), we have also 

{g, V • V’^P 2 )l 2 (o) = -(Vg, Vp 2 )l'^(q) < -(Vg, ^P 2 )l'^{q.) + e||g||_f/i(n) 

= (5) V • Vp2)L2(n) -I- e||g||fl-i(n) < ||5ll/3L2(n)lbi||/3-iff2(Q) -P e|b|l_f/i(o) 
and, by combining these estimates, 


b, V • ypj^)L^(n) < ||5llHi(n)+/3L2(n)|bjr||//i(n)+/3-ifl-2(Q) -p e|b|l/i-i(n) 


By combining these results, we immediately obtain (8.7). This finishes the proof. □ 
Proof of inequality (|5.4[). First note that 


W’f’kPkWn^iil) ~ ||V(bfcPfe)||L 2 (Q) -p \\tpkPk\\L^(n) 

^ ll(VV'fe)Pfc||L2(n) + ||bfeVpfc||L2(Q) -p WipkPkh^ia) 

< IIVbfe||L~(n)lbfc||L2(n) + ||V’fellL°°(n)l|Vp/c||L2(n) + llbfclU“(n)lb/c||L2(n) 

holds due to the product rule and Holder’s inequality. Here, |j • ||loo(o) is the standard 
L°°-norm (essential supremum). As tpkiO is piecewise linear and bfc(C) ^ [Ob]i 
|Vbfe(C)l ^ ^k^ holds for almost all ^ G ft. Using this, |bfc(C)l ^ 1 a standard 
inverse inequality for estimating |lVpfc|b 2 (Q) from above, one obtains HbfcPfellffqn) ^ 
(1 + h^^)|bfc||L 2 (n) and, assuming hk < 1, the desired result. □ 

Proof of inequality (5.5). Assume that the mesh consists of the elements % for 


i = 1,..., Mk- Note that it suffices to show ( |5.5[ ) for each element Ti or that 

{Pk,Pk)L^{Ti) ^ C^{Pk,'f’kPk)L2{Ti) (8-8) 

for all pk G Pk- Note that = 1 on all interior elements (no vertex is located on 


dfl), so (8.8) holds with C = 1 in this case. So, it remains to show the inequality for 
the elements elements where one (case 1) or two vertices (case 2) are located on dil. 
(Note that by assumption each element has at least one vertex in the interior of U). 

Note that, both pk and 'ifk, are linear functions. The integrals in (8.8) can be 
computed using the reference element A := {(^ 1 ,^ 2 ) G (0) 1)^ ^ Ci + C 2 < !}• 

Using the linear transformation : A ^ the substitution rule states that ( |8.8[ ) is 
equivalent to 


|detV$i|(pfc,pfc)i2(7-.) < C'|detV4>i|(pfc,b/cPfc)L2(7-.), 


(8.9) 


where p{^) = u{^~^{^)) and b(0 = Here, |detV$i| is the absolute value 

of the Jacobi determinant of the transformation. As this is a positive constant, it can 
be canceled out. 

For the case 1, the function tpk takes the value 1 on two vertices and the value 0 
on one vertex. We assume without loss of generality, that it takes the value 0 on (0,1). 
So, we obtain ipkiO = 1 — ^ 2 - Using the ansatz Pk{£,iA 2 ) ■= o-o + + 02 ^ 2 , we can 

compute both integrals in (8.9) and obtain that (8.9) holds for C = | (6 + -s/fi). 

For the case 2, the function ifk takes the value 0 on two vertices and the value 
1 on one vertex. We assume without loss of generality that it takes the value 1 on 
the vertex (0,1). So, we obtain bfc(C) = ^ 2 - Here we obtain that (8.9) holds for 
C = (4 -P -y/fi) • This finishes the proof. □ 
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